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ABSTRACT 



Context. The detection of y-rays from dark matter (DM) annihilation is among the scientific goals of the Fermi Large Area Telescope 
(formerly known as GLAST) and Cherenkov telescopes. 

Aims. In this paper we investigate the existence of realistic chances of such a discovery selecting some nearby dwarf spheroidal galax- 
ies (dSph) as a target, and adopting the DM density profiles derived from both astronomical observations and A'-body simulations. 
We also make use of recent highlights about the presence of black holes and of a population of sub-subhalos inside the Local Group 
(LG) dwarfs for boost factor studies. 

Methods. We study the detectability with the Fermi-LAT of the y-ray flux from DM annihilation in four of the nearest and highly 
DM-dominated dSph galaxies of the LG, namely Draco, Ursa Minor, Carina, and Sextans, for which the state-of-art DM density 
profiles were available. We assume the DM is made of Weakly Interacting Massive Particles such as the Lightest Supersymmetric 
Particle (LSP) and compute the expected y-ray flux for optimistic choices of the unknown underlying particle physics parameters. We 
then compute the boost factors due to the presence of DM clumps and of a central supermassive black hole. Finally, we compare our 
predictions with the Fermi-LAT sensitivity maps. 

Results. We find that the dSph galaxies shine above the Galactic smooth halo: e.g., the Galactic halo is brighter than the Draco dSph 
only for angles smaller than 2.3 degrees above the Galactic Center We also find that the presence of a cusp or a constant density 
core in the DM mass density profile does not produce any relevant effect in the y-ray flux due to the fortunate combination of the 
geometrical acceptance of the Fermi-LAT detector and the distance of the galaxies and that no significant enhancement is given by 
the presence of a central black hole or a population of sub-subhalos. 

Conclusions. We conclude that, even for the most optimistic scenario of particle physics, the y-ray flux from DM annihilation in the 
dSph galaxies of the LG would be too low to be detected with the Fermi-LAT. 
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1. Introduction 

Since the first evidence of the presence of dark matter (DM) 
in the universe, scientists have worked to understand its nature 
and distribution. This investigation involves different fields of 
research such as particle physics, cosmology and, observational 
astronomy (e.g. |Bahcall et al. (19 99), Sper gel et al. (2003)[ ). 

The Fermi Large Area Telecope (Fermi-LAT) will 
test theories in which DM candidates are the Lightest 
Supersymmetric Particles (LSPs) such as the neutralinos, 
arising in Supersymmetric extensions of the Standard Model of 
particle physics (SUSY), or the Lightest Kaluza-Klein Particles 
(LKKPs) such as the B^^h, first excitation of the hypercharge 
gauge boson in theories with Universal Extra Dimensions 
(see Bergstrom (2000) [Bertone et al. (20()5)| and references 
therein). Typical values for the mass of these candidates range 
from about 50 GeV up to several TeV. 

Cosmological models, mainly based on A^-body simulations 
in a A-Cold Dark Matter (CDM) framework, successfully re- 
produce relevant characteristics of the universe such as the cos- 
mic microwave background anisotropy and the large scale struc- 



ture of the universe. They also predict well-defined properties 
of DM haloes, whose radial mass density distribution follows 
a universal law, and it is described by a steep power law for 
a wide range of masses ranging from dwarf galaxies to galaxy 
clusters (see, e.g., [Navarro et al. (199 6), Nava rro et al. (1997) 
)[ JMoore et aTTT998), ,Moore et aL (1999) 



Navarro et al. (200' 



Diemand et al. (2005)| l. However, the astronomical community 
is still debating whether DM haloes are characterized by a cen- 
tral density cusp. In fact, haloes with a constant density core are 
in most cases preferred to account for the observed kinematics 
of galaxies (see Binney (2004) for a review). 
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Generally speaking, the uncertainty in the choice of the den- 
sity profile can result in several order of magnitude of indetermi- 
nation on the y-ray flux prediction, which already suffers from 
the high uncertainties arising by the unknown underlying parti- 
cle physics (Forne ngo et al. (2004)] l. For this reason it would be 
important to derive the DM density profile of galaxies directly 
from the available kinematic data. Although data-sets for the 
very inner part of the galaxies are scarce and affected by large 
errors, the situation is not better in A^-body simulations, whose 
resolution goes down to 0.05 times the virial radius at most. 
Using real data we have the advantage of deriving a flux pre- 
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diction which takes into account the pecuharity of each galaxy, 
without any model-dependent generalization which would arise 
the astrophysical uncertainties. 

The expected y-ray flux at the telescope from a given source 
is directly proportional to the DM density squared along the line- 
of-sight (LOS), and inversely proportional to the square of its 
distance. The best targets are therefore nearby dense objects such 
as the local dwarf spheroidal galaxies (see Mateo (1998)| for a re- 
view). Indeed, in the last decade, the large collecting area of the 
8-m class telescope and the use of multi-fiber spectrographs al- 
lowed astronomers to obtain high-resolution spectra of a large 
number of stars. This made possible to isolate the galaxy mem- 
ber stars, to measure their radial velocity with an accuracy of few 
km s ' and to build accurate dynamical models of a number of 
such systems (see| Gilmore et al. (2007 )' and references therein). 

Annihilation of y-rays in dSph galaxies would give a 
clean signal because of the absence of high astrophysi- 
cal uncertainties in modeling the expected background and 
could hopefully be detected with upcoming experiments 
like the Fermi-LAT. Many authors studied the feasibil- 
ity of such a detection, using a large variety of cuspy 
and cored universal density profiles, reflecting the theo- 
retical as well as the experimental uncertainties. Diff'erent 



works ( 


Baltz et al. (2000) 


Tyler (2002) Peirani et al. (2004)| 


IPieri & 


Branchini (2004) 


Bergstrom et al. (2006) I found that 



only the presence of a spike and/or an enhancement due to 
dumpiness and/or a more favourable combination of the un- 
known particle physics parameters could make the Draco dwarf 
galaxy observable with the Fermi-LAT. Strigari et al. (2007)| are 
optimistic about the detection of Draco with the Fermi-LAT 
in 5 years. They adopted a King profile (King (1966)] l to de- 
rive the surface density of the stellar luminosity. This was de- 
projected and converted into the stellar mass density by adopt- 
ing the typical range of the mass-to-light ratio of dSphs. The 
luminous mass they derived is at the very least one order of 
magnitude below the mass of the DM halo. For the halo they 
assumed a NFW mass density profile ( [Navarro et al. (1996)[ 
[Navarro et al. (1997)| l whose free degenerate parameters were 
constrained by marginalizing over the stellar velocity disper- 
sion anisotropy parameter. |Colafrancesco et al. (2007)[ showed 
how diffiise radio emission would actually be a more promis- 
ing process to look at in order to detect a DM signal. They also 
claimed that he presence of a supermassive black hole (SBH) at 
the centre of Draco, which could enhance the y-ray signal up to 
detection, is not actually excluded by experiments. Detection of 
annihilation y-rays from Draco has been recently excluded by 
[Sanchez-Conde et al. (2007)1 through the use of density profiles 
which are compatible with the latest observations. 

In this paper we use the latest available astrophysical 
measurements for four of the nearest and highly DM-dominated 
dSph galaxies of the Local Group, namely Draco, Ursa Minor, 
Carina, and Sextans, to compute the expected y-ray flux from 
DM annihilation. 

In Sect. [2] the most optimistic particle physics scenarios 
and the DM density profiles derived both from the available 
kinematic measurements and from A^-body simulations are 
used to predict the expected y-ray flux from DM annihilation 
in Draco, Ursa Minor, Carina, and Sextans. In Sect. [3] the 
predicted flux is compared with the experimental sensitivity 
of the Fermi-LAT. The presence of DM clumps and a central 
SBH could enhance the y-ray flux. But their effects have to 
be rescaled for the limits imposed on the extragalactic y-ray 
background (EGB) by the Energetic Gamma-Ray Experiment 



(EGRET) and on the y-ray flux in Draco by the measurements 
of the Major Atmospheric Gamma-Ray Imaging Cherenkov 
(MAGIC) telescope. Our conclusions are given in Sect. [4] 

The main diff'erences with the other papers that already dis- 
cussed this argument are the following: we show that, even 
adopting a very lucky case for the unknown particle physics sec- 
tor, the expected flux from the DM halo is about two orders of 
magnitude below the detectability limit of the Fermi-LAT exper- 
iment; we also show how the use of a cored or a cuspy profiles 
does not produce any relevant eff'ect in the expected y-ray flux, 
because of a combination of the galaxy distance and the angular 
acceptance of the Fermi-LAT; we then show that the current lim- 
its on the mass of the supermassive black hole (SMBH) inside 
Draco lead to an insignificant boost factor due to the presence 
of such a SMBH; and we numerically compute the boost fac- 
tor due to the presence of a population of subhaloes inside the 
dwarf galaxies, limiting the possible range of models for the sub- 
subhalo structure making use of the constraints imposed by the 
EGRET extragalactic measurements; the boost factor due to the 
presence of sub-subhaloes is computed in two ways: first, we ob- 
tain it integrating over the whole volume of the galaxy, as done, 
e.g., in |Strigari et al. (2007)| This gives the correct boost factor 
when considering cosmological haloes. But, if we consider the 
closer dwarf galaxies, we have to take into account that only the 
very inner part of the galaxy is observed within the angular reso- 
lution of the instrument. We therefore also compute the angular 
dependence of the boost factor due to sub-subhaloes. Although 
we find a huge enhancement of the flux far from the galaxy cen- 
ter, there is actually no enhancement along the LOS pointing 
toward the galaxy center. As a last improvement with respect to 
the other papers, we compare our predictions with an recently 
released detectability map for the Fermi-LAT which takes into 
account the response of the detector to different energies and in- 
cidence angles, as far as effective energy and angular resolution 
are concerned. 



2. y-ray flux from dark matter annihilation 

The y-ray flux <t)y from DM annihilation can be factorized into a 
term Opp involving the particle physics and a term Ocosmo where 
astrophysics, cosmology, and experimental geometry play the 



main role. It is 

%{Ey, l/r, AO) = <^Pp(Ej.) X Oeosmo(l/'. AO). 

The particle physics factor is given by 
Opp(£,).l^yB,. r 



(1) 



(2) 



where otdm is the DM particle mass, cr3„n is the annihilation 
cross section, and v is the relative velocity. cr.^n„v determines 
the number of annihilations. Bj is the branching ratio into a fi- 
nal state /. It represents the probability that the final state / is 
the result of one annihilation. dN^/dEy is the yield of photons 
produced by the final state / in one annihilation, and Eth is the 
threshold energy above which the flux is computed. So far no 
assumptions have been made on the nature of the DM particles. 
For a complete set of the allowed values for the previous param- 
eters the reader is referred to Fornengo et al. (2004) 
The astrophysical/cosmological factor is given by 



□((A, AO) = f f d<f>de I dA 
J Jaq Jlos 



Pdm('-) 



^2 



(3) 
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where t]/ is the angle of view from the halo centre which de- 
fines the LOS and AQ corresponds to the angular resolution 
of the instrument. It is in principle function of the photon en- 
ergy, Ey, though in the following we will assume for simplicity 
AQ = 10"^ sr (corresponding to a cone of view with angular 
opening of 0.1 degree along the LOS). This corresponds to the 
geometrical acceptance of the Fermi-LAT detector. Indeed, we 
made the simplifying assumption that the angular resolution of 
the Fermi-LAT is 0.1 degree over the entire energy range. This 
resolution is reached only for 10 GeV photons and for incidence 
angles less than 50°. A is the coordinate along the LOS, and r the 
radial coordinate inside the halo, pouif) is the DM mass density 
profile, which is a factor of primary importance when deriving 
the y-ray flux from DM annihilation in a given source. Finally, 
J{x, y, z\A, 9, <p) is the Jacobian determinant. 

2. 1 . The particle physics factor 

In Fig. [U we draw the factor Opp integrated above a threshold 
energy Zith as a function of Eth. Shown in the plot is the re- 
sult of the computation for a 40 GeV, 100 GeV and 1 TeV DM 
particle annihilating into quarks bb. For the photon yields we 



have used the parametric formula from Fornengo et al. (2004) 



and introduced the pion bump feature at low energies. For 
each mass, we adopted the most optimistic value for iXannV, as 
it is computed with DARKSUSY ( Gondolo et al. (2004)| l and 
allowed by WMAP-nSDSS measurements (e.g, see Fig. 3 in 
|Pieri et al. (2008)| ). However, really few models lie in those very 
fortunate parts of the phase-space. In details, we used o-.^„„v - 



X 10 



* cm-* s^'for a 40 GeV DM particle, 0"^, 
s^'for the 100 GeV particle case and cr^„, 
-If 



= 10 
= 10 



-25 
-26 



We do 



not 



consider here the result of 



Bringmann et al. (2007) 



who pointed out how so far ig- 
nored efl'ects of electromagnetic radiative corrections to all 
leading annihilation processes in the Minimal Supersymmetric 
Model or in the Minimal SUperGRAvity mediated supersym- 
metry breaking scenarios can induce a y-ray flux enhancement 
up to three-four order of magnitudes with respect to the y-ray 
secondary flux produced in the annihilation cascade, when 
integrating over energies greater than 60 % of otdm, even for 
LSP masses well below the TeV scale. A careful study of 
the eff^ect of internal bremsstrahlung would be interesting for 
instruments with higher sensitivity at higher energies, such as 
Cherenkov telescopes, and is beyond the goal of this paper 

In the following, we will refer to a 40 GeV DM particle with 
o'annV = 3 X 10~^^ cm^ s~\ annihilating into bb as to our best 
case scenario when studying the maximal y-ray flux prediction 
integrated above 100 MeV The reader must be aware that this is 
not the most likely model, and that the real flux could be orders 
of magnitude smaller. 



2.2. The astrophysicai/cosmoiogical factor 

In this section we derive the value of <I>cosmo for four dSph galax- 
ies of the Local Group both from the state-of-art DM density 
profiles available in literature and from CDM A^-body simu- 
lations. Their positions, masses, and distances are reported in 
Table [U 



- Draco: [Gilmore et al. (2007)] calculated the DM density ra- 
dial profile of the Draco dwarf galaxy. It was derived by 
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Fig. 1. Behaviour of ^pp{> E^h) as a function of E^h for different 
models of the DM particle, computed for a: 40 GeV (solid line), 
100 GeV (dotted) and 1 TeV (dashed) DM particle annihilating 
into bb. The value of cra„nV' has been chosen as representative of 
the best value for that mass, as it is computed with DARKSUSY 
and allowed by WMAP-nSDSS measurements (see details in the 
text). 

Table 1. The sample galaxies. Distances and masses are taken 
from |Gilmore et al. (2007)1 



Object 


/ 


b 


Mass 


Distance 




deg 


deg 


Mo 


kpc 


Draco 


86.37 


34.77 


2.8 X 10' 


80 


Ursa Minor 


104.98 


44.86 


1.3 X 10^ 


66 


Carina 


260.09 


-22.28 


1.3 X 10' 


87 


Sextants 


243.42 


42.16 


1.9 X 10' 


80 



[Wilkinson et al. (2004)| from the radial profiles of the veloc- 
ity dispersion and surface brightness by solving the Jeans 
equations under the assumption of isotropic orbital structure. 
The velocity dispersion radial profile extends out to about 35 
arcmin from the centre (corresponding to 0.8 kpc). It is char- 
acterized by an almost constant value of about 13 km s"', 
with a decrease to about 5 km s"' at the last observed radius. 
The available data allowed to derive the mass density profile 
of the DM between about 0.1 and 0.5 kpc from the galaxy 
centre. The mass density increases out to the innermost ob- 
served point (Fig. |2]l. 

Recently, an independent mass density profile for Draco has 
been obtained bv lPefiarrubia et al. 20081 They used the data 
by lWilkinson et al. ( 2004)'and pufioz et al. (2006 )lto recon- 
struct the mass distribution of the galaxy. They assumed that 
the galaxy is composed by a luminous component described 



by means of a King model King (1966), and a DM compo- 



nent described by a NFW model. In this way they derived 
the concentration parameter of the DM halo component di- 
rectly from a fit to the data, instead of assuming it from the 
CDM cosmology. A total mass of 6.2 x 10'' M© was found, 
which is somehow larger than expected for dSph galaxies 
( |Mateo (1 998)1 . This DM density profile is shown in Fig.|2] 
The same procedure was applied to the other galaxies we 
analyze. Yet, as we will show in the following, this profile 
implies a more pessimistic y-ray flux prediction. Since we 
are interested in the most optimistic scenarios which possi- 
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Draco 



The radial 
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-0.6 



-0.3 



profile 

by " 



Gilmore et al. (2007) 



mass density in 
(solid line), 



ILokas et al. (2005)| (dotted ), W alker et al. (2007) | (long-dashed), 
and Pen arrubia et al. 20081 (dashed). Also shown are the den- 
sity profiles derived from numerical simulations, namely the 
standard NFW (long-long-dashed) and the Einasto (long-dot- 
dashed) radial profiles. 



bly could lead to detection, we will consider the DM profiles 
derived by Penarrubia e t al. 20081 only in the case of Draco 
to derive the model uncertainty, while we will no further con- 
sider it for the other galaxies. 

A third DM profile for Draco was obtained by 
They used the data-set by 
and assumed a Sersic law 
the distribution of the lumi- 



Lokas et al. (2005) 



Wilkinson et al. (2004)| 



( [Sersic (1968)| ) to describe 
nous component. Concerning the DM density distribution, 
they assumed a modified NFW with an inner cusp and an 
exponential cut-off to take into account a possible tidal strip- 
ping in the outer regions of the galaxy. A tidal interaction 
does not affect the DM mass density profile in the centre, but 
produces a mass loss for radii larger than the so-called break 
radius ( [Kazantzidis et al. (2004)l i. [Lokas et al. (2005)] break 
the degeneracy between the mass distribution and velocity 
anisotropy by fitting both the LOS velocity dispersion and 
kurtosis profiles. They found a total mass of 7 x 10^ Mq. 
The coiTesponding radial profile of the DM mass density 
profile is shown in Fig.[2| 

Finally, we show in Fig. [2[the density profile obtained by 
[Walker et al. (200T) adopting a one-component King profile 
and a NFW profile with constant anisotropy parameter for 
the luminous and DM components, respectively. 
A cored profile seems to be preferred for the DM mass 
density when no parametric function is imposed in fitting 
the data. A primordial density core would exclude a pure 
CDM scenario, rather pointing toward a warm dark matter 
particle. Yet, there are different studies about the possi- 
bility of dynamically remove the CDM cusp in the dwarf 
galaxies, involving phenomena such as stellar feedback 
(e.g., [Mashchenko et al. (2006)[ [Read & Gilmore (2005)| l 
or dynamical friction of DM/baryons subhaloes 
( [Romano-Dfaz et al. (2008)) . The reader should be 
aware that the topic is still controverse (see, e.g, 
[Gnedin & Zhao (2002)[ l and there is no univoque con- 
sensus about the realistic possibility that CDM cusps in 
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Fig. 3. The radial profile of the DM mass density in Ursa Minor 
(long-dashed line), Draco (solid line). Carina (dotted line), and 
Sextans (short-dashed line) as derived by [Gilmore et al. (2007)[ 



dwarfs may be reduced to a core. 

In lack of negative evidences we keep on using cored 
profiles associated with CDM particles in our discussion. 
Although the main aim of the present paper is to present 
results based on density profiles directly inferred by 
astronomical data, it is worth superimposing on Fig. [2| 
the density profiles derived from numerical simulations. 
[Stadel et al. (2008)] have recently obtained from A^-body 
simulations a best fit to a MW-sized halo which is a simple 
power law in dlog(p)/dlog{r) (called the Stadel & Moore 
profile -S&M-). In the lack of halo mass scaling relations 
for the parameters of the S&M profile, we show only the 
NFW and the Einasto profile computed for a 10^ Mq halo 
at a distance of 80 kpc, with the concentration parameter 
given by Kuhlen et al. (2008)[ (c = 20.2, r, = L02 kpc, 
pf'''^ = S^exWM^^c^f™"" = 8.9x lO^Mokpc^l) 
Ursa Minor: The DM mass density profile for the Ursa 
Minor dwarf galaxy was taken from Gilmore et al. (2007)[ It 
was derived by Wilkins on et al. (2004) ' from the radial pro- 
files of the velocity dispersion and surface brightness by 
solving the Jeans equations under the assumption of isotropic 
orbital structure. The velocity dispersion radial profile ex- 
tends out to 45 arcmin from the centre (corresponding to 
0.9 kpc). It is characterized by constant value of about 12 
km s"', showing a sharp drop to about 2 km s"' only at the 
farthest observed radius. The data allowed to derive the DM 
density distribution in the radial range between about 0. 1 and 
0.5 kpc. It is similar to that of Draco (Fig.[3]l. 
Ursa Minor w as also studied in [Strigari et al. (20 07) using 
the data by jPalma et al. (2003)[ The light distribution 
was derived considering a two-component, spherically- 
symmetric King profile. They used a Jeans equations and 
adopted a NFW DM halo to the derive the radial profile 
of the velocity dispersion fitting the data. The anisotropy 
parameter /3 has been empirically set to the value of 0.6. 
In this paper we consider the two NFW models of the DM 
density profile given by [Strigari et al. (2007)] Model A 

Model B has 



has = 0.63 kpc and ps = 10^ M© kpc 
Ks - 3.1 kpc and pj = 10^ Mokpc"^. 
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Carina: |Gilmore et al. (2007)| calculated the DM density 
radial profile of the Carina dwarf galaxy. It was derived 
from the radial profiles of the velocity dispersion and 
surface brightness by solving the Jeans equations under the 
assumption of isotropic orbital structure. The available mea- 
surements extend out to about the tidal radius of the galaxy, 
which corresponds to about 25 arcmin (corresponding to 0.6 
kpc). The velocity dispersion is characterized by a constant 
value of about 8 km s The DM mass density profile is 
derived out to 60 pc from the centre and it shows a constant 
density core (Fig.[3]l. 

Sextans: The DM mass density profile for the Sextans dwarf 
galaxy was taken from Gilmore et al. (2007)] It was de- 
rived by [Wilkinson et al. (2006)| from the radial profiles of 
the velocity dispersion and surface brightness measured by 
Kleyna et al. (2004) by solving the Jeans equations under the 
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assumption of isotropic orbital structure. The velocity dis- 
persion radial profile extends out to about 47 arcmin from 
the centre (corresponding to 1.1 kpc). It is characterized by 
a constant value of about 8 km s"', with a possible decrease 
to about 3 km s"' at the last observed radius. The available 
data allowed to derive the mass density profile of the DM be- 
tween about 0.2 and 0.8 kpc from the galaxy centre. It shows 
a constant density core (Fig.O. 

We would like to underline that the experimental results 
which we will use in our analysis give both cuspy and cored 
profiles. Although CDM simulations predict only cuspy haloes, 
we will keep on considering cored profiles and CDM particle be- 
cause there may be mechanisms of gravitational heating of the 
dark matter by baryonic components which could reconcile ob- 
servations of cored density profiles with the central density cusps 
of the CDM predictions. 

Eq. [3] has then been integrated along the LOS adopting the 
DM density profiles we derived for each dSph galaxy. The re- 
sult of this integration for the four different profiles inferred by 
the data as well as for the two profiles derived from numeri- 
cal simulations for the Draco dSph galaxy is found in Fig. |4] 
The behaviour of these curves reflects the different DM distri- 
butions shown in Fig. |2l The fit to the data was obtained in 
the radial interval between 80 and 630 pc. The DM mass den- 
sity profiles are extrapolated in the innermost and outermost 
galaxy regions. At large radii the DM mass density derived by 
IPeflarrubia et al. 20081 who adopted a NFW density profile with 
no tidal disruption, is higher than that by [Lokas et al. (2005)| 
and [Gilmore et al. (2007)1 Indeed, it is comparable to the re- 
sult obtained by 'Walk er et al. (2007)| who also fit a NFW pro- 
file. Actually, the total mass derived by IPeflarrubia et al. 20081 
and [Walker et al. (2007)] is larger than the one typically found 
for this kind of galaxies ( [Mateo (1998)[ ). This behaviour reflects 
in the radial trend of the corresponding <I>cosmo (Fig.|4|, which is 
higher at large radii with respect to those based on the results by 
[Lokas et al. (2005)[ and [Gilmore et al. (2007)[ The results based 
on the [Walker et al. (2007)[ profile give a higher value of Ocosmo 
in the inner galaxy with respect to the 'Peflarrubi a et al. 20081 
since the former predict a larger mass content at small radii (see 
Fig. 13. 

The low DM mass density observed at large radii in the 
NFW profile of [Lokas et al. (2005)[ is due to the mass stripping 
induced by a tidal interaction. Their DM density is higher in the 
centre, while the cored density profile by Gilmore et al. (2007)] 
allocates more mass at intermediate radii (Fig. [2|). Though 
biased by the different derived masses, this effect is due to mass 
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Fig. 4. The astrophysical/cosmological contribution <l'cosmo to 
the y-ray flux derived as a function of the angular dis- 
tance from the galaxy centre for Draco from the DM mass 
density radial profiles by Gilmore et al. (20()7)| (solid line). 



[Lokas et al. (2005)| (dotted ) ,Walker et al. (2007)[ (long-dashed), 
and IPeflarrubia et al. 20081 (dashed). Also shown are the results 
obtained from the density profiles derived from numerical sim- 
ulations, namely the NFW standard (long-long-dashed) and the 
Einasto (long-dot-dashed) profiles. 



conservation since the two models have about the same tidal 
radius. The DM mass density profile by Gilmore et al. (2007)1 
gives a larger <I>cosmo for radii larger than 0.1 degree. At smaller 
radii it gives the same contribution as the DM mass density 
profile by [Lokas et al. (2005)[ (see Fig. lU. The Einasto profile, 
which predicts more mass at the intermediate radii which are 
resolved by the angular resolution of 0.1 degrees, gives the 
highest value of Ocosmo, while the NFW profile gives the same 
contribution as IPeflarrubia et al. 20081 The values of the results 
for the NFW and Einasto profiles depend on the mass adopted 
for the computation. We have used here lO'^M© because the 
relative density profile was compatible with the amplitude of the 
profile inferred by the data. 

It is worth stressing here that the mass modeling of Draco 
produces only a difference of a factor 2 to 3 in the flux predic- 
tions, while the indetermination arising from the unknown parti- 
cle physics can sum up to several orders of magnitude. 

To investigate the reason why, e.g., the cuspy pro- 
file by jLokas et al. (2005)[ and the cored profile by 



[Gilmore et al. (2007)[ give the same value of Ocosmo to- 



wards the centre of Draco, we considered a Draco-like dSph 
galaxy and changed its distance from the observer. We then 
computed Ocosmo toward the centre of the galaxy. The result for 
the two profiles is plotted in Fig.|5las a function of the imposed 
distance. The closer is the galaxy, the greater is the contribution 
to Ocosmo due to the cuspy radial profile of the DM mass density. 
The geometrical acceptance of the Fermi-LAT detector is able 
to resolve the central cusp of the galaxy only if this is located at 
distances smaller that 90 kpc. Further out, the two profiles give 
more or less the same result. Curiously enough, the true location 
of the Draco dSph (80 kpc from us) lies exactly at the border 
of this region, so that we can conclude that no matter whether 
we choose either the cuspy DM profile by [Lokas et al. (2005)1 
or the cored profile by [Gilmore et al. (2007)[ the estimate of the 
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Fig. 5. The astrophysical/cosmological contribution Ocosmo to 
the y-ray flux computed in the centre of a Draco-like galaxy as 
a function of its distance to observer. The DM was considered 
to be radially distributed as in |Gilmore et al. (2007)| (solid line) 
and [Lokas et al."(2005), (dotted) 
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Fig. 6. The astrophysical/cosmological contribution Ocosmo to 
the y-ray flux derived as a function of the angular distance 
from the galaxy centre for Ursa Minor from the DM mass den- 
sity radial profiles by 'Gil more et al. (2007)| (long-dashed line), 
IPenarrubia et al. 20 08 (dashed), and for the two fit to the NFW 




profile proposed in Strigari et al. (2007) (solid and dotted) 



amount of y-rays expected from DM annihilation in the central 
region of the galaxy will not change. 

In Fig. |6] we plot the value of Ocosmo for Ursa Minor for 
the cored iGilmore etal. (2007)| profile, as well as for the cuspy 
IPenarrubia et al 20081 and for th e two fit to flie NFW profile 
proposed in 'Strigari et al. (2007)1 As in the case of Draco, the 
[Peiiarrubia et al. 2008 profile gives the lowest value, while the 
two NFW models of [Strigari et al. (2007)| bracket the cored 
value at small angles. 

What has been afore-discussed holds for the other dSph 
galaxies considered in this analysis too. As an example, 
in Fig. [T] we plot the value of Ocosmo obtained using the 
[Gilmore et al. (2007)[ profile for the four dSph galaxies consid- 
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Fig. 7. The astrophysical/cosmological contribution Ocosmo to 
the y-ray flux derived as a function of the angular distance 
from the galaxy centre for Ursa Minor (long-dashed line), 
Draco (solid line). Carina (dotted line), and Sextans (short- 
dashed line) from the DM mass density radial profiles by 
[Gilmore et al. (2007)[ 



ered in this analysis. The values obtained using cuspy profiles 
will not deviate significantly from these values. 

Finally, FiglSj shows the values of Ocosmo obtained the pro- 
files by Gilmore et al. (2007) 1 and computed for the LOS point- 
ing toward the centre of the four dwarfs. These values are com- 
pared to the curve obtained for the smooth halo of the MW, 
obtained using Eq|3] an angular resolution of 0. 1 degrees and 
the NFW profile for the MW (Mmw = IO'^Mq, c = 7.55, 
r, = 27.3 kpc). We observe that the dSph galaxies shine above 
the smooth Galactic halo at their position in the sky. Even more, 
we can say that, e.g. Draco is brighter than the Galactic halo at 
all angles larger than 2.3 degrees above the Galactic center 
The central values of Ocosmo for the Sagittarius dwarf galaxy and 
Large Magellanic Cloud (LMC) are shown in Figl8]for a sake 
of comparison with those of the other dSph galaxies we stud- 
ied in detail. The Sagittarius dwarf galaxy is located at a dis- 
tance of about 24 kpc. Although it is heavily interacting with 
the MiUcy Way, it has a surviving stellar component thus it is 
likely to have a surviving dark matter halo. The observations 
suggest that it is dark matter dominated with a central stel- 
lar velocity dispersion of about 10 km s ' [Ibata etal.(1997)| 
According to recent observations and semi-analytic modelling 
(e.g. [Strigari et al. (2008)l[Macci6 et al. (2008)) , the data consis- 
tent with all the DM halo of the dSph galaxies lie in the range 
between 20 and 40 km s We then modeled the inner regions 
of the DM halo of the Sagittarius dwarf with the same scale pa- 
rameters as Draco (see |E vans et al. (2004)) l by assuming a NFW 
mass density profile and a mass of M = 10^ Mq. The LMC is 
located at about 50 kpc. We adopted for its DM halo the stripped 
NFW profile used by|Tasitsiomi et al. (2004),(M ~ 10"' Mq). 



3. Predictions for observation with tlie Fermi-LAT 

The map of the Fermi-LAT sensitivity to point sources of DM 
annihilations has been obtained by [Baltz et al. ( 2008^ using the 
released Fermi-LAT response functions. The sensitivity map was 
obtained for 55 days of observation and it shows the minimum 
flux above 100 MeV which is necessary in order to achieve a 5cr 
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Fig. 8. The astrophysical/cosmological contribution Ocosmo to 
the y-ray flux derived as a function of the angular distance from 
the center of the Milky Way centre, computed for a MW NFW 
smooth halo (solid curve), and for the central angular bin of Ursa 
Minor (filled circle), Draco (open circle). Carina (open triangle), 
and Sextans (filled triangle) derived using the DM mass density 
radial profiles by |Gilmore et al. (2007)| Also superimposed are 
the values for a NFW fit to the Sagittarius and the LMC galaxies. 



detection. The significance of the observed signal given the local 
background counts is assigned with a maximum likelihood anal- 
ysis assuming Poisson statistics. |Baltz et al. (2 008) found that 
the sensitivity suffers very little dependence from the underlying 
particle physics. Therefore, the obtained values can be consid- 
ered valid as long as the source appears point like in the sky, that 
is as long as its angular size does not exceed 0.25 degrees. As 
it can be seen in Fig. |7] the y-ray flux expected from our dSph 
galaxies decreases by almost one order of magnitude at the an- 
gular distance of 0.25 degrees from the galaxy centre. For this 
reason, we can assume they are point-like sources and use the 
results of |Baltz et al. (2008)| for reference. Draco, Ursa Minor, 
and Sextans lie in a region of the sky where the 5cr detection 
flux above 100 MeV is 1.5 x 10"^ ph cm"^ s^'in ~ 2 months. 
This translates into 0^^^^"^ = 6 x 10"' ph cm"^ s"' in 1 year of data 
taking and in the units we used throughout this paper. In the case 



■ lyr 



= 8 X 10 ph cm ^ s ^ since the galaxy is 



of Carina, it is y,^ 
closer to the Galactic plane. 

If we consider the best value for <I>pp (> 100 MeV) from 
Fig.[I]((tpp ~ 6x lO-** cm'* kpc"' GeV"^ s"' sr"') and the 
average value of <l>cosmo toward the galaxy centre (tf/ - 0) from 
Fig. |4] and Fig. |6] we end up with the following best-particle- 
physics-case estimates for the y-ray flux from DM annihilation 
in Draco: 



Draco / 



y (> 100MeV) = (4.6 + l.l)xlO""phcm"^s" 
and Ursa Minor: 



(4) 
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Fig. 9. Differential y-ray fluxes as a function of the energy, ex- 
pected from DM annihilation in the center of Ursa Minor. Fluxes 
are computed using the best value for Ocosmo given by model A 
of Strigari et al. (2007) , and are presented for a 40 GeV (solid 
line), 100 GeV (dotted) and 1 TeV (dashed) DM particle annihi- 
lating into bb. 



The error is given by the standard deviation for the values of 
Ocosmo (lA - 0) obtained using different DM density profiles in- 
ferred by dynamical modeling. We therefore do not focus on a 
specific profile when giving the value of the expected y-ray flux. 
Indeed, our result is obtained by averaging over different fits to 
the data. It is worth noticing that predictions made using profiles 
inferred by astronomical data are robust within a 60 % relative 
error which is expected while changing fit. This means that they 
can provide a reliable order-of-magnitude estimate of the real 
flux. We will not further consider the case of Carina and Sextans 
since they give a lower flux. Yet, the calculations of the expected 
y-ray fluxes from DM annihilation in these galaxies are straight- 
forward. 

Even in the case of Draco and Ursa Minor the upper value 
of the predicted flux within the error ends up to be 2 orders of 
magnitude below that required for detection in 1 year of data 
taking with the Fermi-LAT, referring to |Baltz et al. (200¥ ) This 
means that there is no hope of detection unless we allow for the 
presence of boost factors. 

Though brigthest than the dSph considered in this analysis, none 
of the Sagittarius dwarf galaxy and LMC have a predicted flux 
which could be detected with the Fermi-LAT. 

In Fig.|9]we show the differential y-ray fluxes expected from 
DM annihilation in the center of Ursa Minor for a 40 GeV, 100 
GeV and 1 TeV DM particle annihilating into bb. Fluxes are 
computed using the best value for Ocosmo given by model A of 

anni' havc been chosen as 



Strigari et al. (2007) The values of cr^, 
in FigUJ 



We note that if we use the same values for the annihilation 
cross-section and for the mass ((r^nnV = 5 x 10"^^ cm^ s"', 
- 46 GeV) as in Strigari et al. (2007)' as well as their model 
A for the density profile, we find a prediction for Ursa Minor 
which is ~ 10 times smaller than their value. In fact, we get 
0^(> 5 GeV) ~ 2.5 [5.4] x lO^'^ phcm^s"' for annihilation 



0™'"(> 100 MeV) = (6.0 ± 3.8) x 10"" phcm"'' s 



1-n 



-2 „-l 



into bb [t^t ] to be compared with their value Oy(> 5 GeV) ~ 
3 X 10"" phcm^s"'. This is due to the over-estimated num- 
(5) ber of photon yields above 5 GeV (^'"^^^ dEdNy/dEy ~ 4.2) 
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which is derived in their paper. We found a number of pho- 
ton yields which is an order of magnitude smaller both us- 
ing the Fomengo et al. (2004) and the Bergstrom et al. (1998) 
parametrization for dNy/dEy, the latter being the one used by 



Strigari et al. (2007) 



Investigating possible sources of astrophysical boost factors 
becomes necessary in order to understand the feasibility of a 
DM signal detection with the Fermi-LAT. To this purpose, in the 
following sections we account for the effect of the presence of 
clumps or of a SBH inside the dSph galaxies we are considering. 

3.1. Boost factor due to the presence of dark matter clumps 

According to the CDM scenario, each halo formed through 
the merging and accretion of smaller haloes, which still sur- 
vive and orbit inside the larger one. The minimum mass of 
these subhaloes is ~ 10"^ M© according to analytical esti- 
mates ( Green et al. (2004)|[Green et al. (2005)| l. High-resolution 
A'-body experiments, though they stop at high redshift (z = 26), 
are able to resolve field haloes as small as ~ 10"^ Mq. Their 
mass function is well approximated by a power law 



dn(M)/dln(M) ocM"", 



(6) 



with a - I, independently of the mass of the host halo, over 
the large redshift range between and 75 and the mass inter- 
val between about 1 0'^ and 10'° Mq ( |Diemand et al. (2005)| 
IGiocoU et al. (2008)| l. 

Assuming that the radial distribution of subhaloes traces that 
of the host galaxy, we can model the number density of sub- 
haloes per unit mass at a distance R from the galaxy centre as 



Psh(M,i?) 



: AM-^0(R 



(7) 



where A is a normalization factor which takes into account the 
hypothesis that 10% of the Milky Way (MW) mass is distributed 
in substructures with masses in the mass range between 10^ and 
10'° Mq dPiemand et 117(20 05) ) The effect of tidal disruption 
is accounted for by the Heaviside step function 0(r - r^in{M)), 
where the tidal radius rmm{M) is computed according to the 
Roche criterion as the minimum distance at which the subhalo 
self-gravity at its scale radius equals the gravity pull of the halo 
host computed at the orbital radius of the subhalo. The debate 
on the survival and partial disruption of these haloes is still 
open, and many issues are still unsolved, such as the true mass 
function after tidal interactions in the host halo, the inner struc- 
ture and concentration of the subhaloes themselves. We refer to 
|Pieri et al. (2008)| for a detailed discussion of the problem. 

Once we assumed a model for the subhalo population, the 
boost factor due to the presence of clumps distributed accord- 
ing to Pf,hiM, r) is computed as the ratio of the integral over the 
galaxy volume of the density squared including subhaloes to the 
same integral for the smooth galaxy only: 



BF: 



4i dVplum + 4i L ^^^^^ Lie dVp^ 



SH 



LdVp: 



2 

galjSm 



(8) 



where Pgai.sm is the smooth profile of the DM component of the 
host halo which is not virialized into clumps and ph is the internal 
DM density profile of each subhalo. 

|Pieri et al. (2008)| found a relationship between the different 
subhalo models leading to more or less impressive boost fac- 
tors for the MW, the total number of photons produced at high 
galactic latitudes by the annihilation of DM particles in all the 
subhaloes falling into a given cone of view (of the order of 10''), 



the EGRET measurement of the extragalactic y-ray background 
(EGB), and the allowed particle physics contribution. They ob- 
served how a given model for the subhalo population can not pre- 
dict a number of photons greater than those observed by EGRET 
at high latitudes, where the y-ray flux is thought to have a dif- 
fuse origin. Consequently, a maximum number of predictable 
photons exists. This means that the two factors Opp and Ocosmo 
must be tuned in order not to exceed the EGRET limit. In the 
most optimistic case, they will be tuned as to give exactly the 
number of photons observed by EGRET. This means that, if we 
assume a subhalo model for the MW, the value of Opp can be 
shifted down or up to match the EGRET level (up to the level of 
the best-particle-physics case of FigdJ. Now, in the lack of accu- 
rate models which account for the presence of subhaloes inside 
subhaloes, we make the simplifying assumption that the subhalo 
population of a dSph galaxy is described by the same subhalo 
model which we have assumed to be valid in the MW, so that the 
restrictions on <t)pp must still hold. 

We have computed the boost factors in the cases of Draco 
and Ursa Minor, for all the subhalo models considered in 
Fie ri et al. (2008)| (PBB08 in FigHOl using Eq.[8] We found that 
the values for the boost factors range from 1.6 to 850, but when 
applying the corresponding limits on Opp, we end up, for any 
clump model, with an estimate of the maximum flux which may 
be produced by the clumps in Draco and Ursa Minor which is 
still compatible with the EGRET EGB and with the constraints 
given by particle physics shown in Fig.[T] The overall maximum 
enhancement of the flux obtained using Eq. [8] once scaled for 
the EGRET limit, is of a factor 70. 

The boost factor due to the presence of subhaloes has 



been computed analytically also in Strigari et al. (2007) and 
[Kuhlen et al. (2008)] The overall values is of the same order 
of magnitude of the one we obtain using the Bzo.ref model of 
Fie ri et al. (2008)| For that specific model, we obtain BFsh - 2, 
which is not expected to give an enhancement of the y-ray flux 
which could be significant for detection. 

The values of the boost factors obtained using Eq. [8] as 
well as those obtained analytically in Strigari et al. (2007)| and 
|Kuhlen et al. (2008)] hold when integrating over all the galaxy. 
They are thus the numbers to consider when the galaxy is so far 
as to be pointlike inside the detector acceptance. This is indeed 
not the case for the nearby dwarfs considered in this analysis. 

To understand what could really be the effect of sub- 
subhaloes in the closest dwarfs, we assumed a NEW profile 
for the substructures, defined by the concentration parameter 
c distributed according to a log-normal probability P{c), and 
have computed their contribution to the annihilation signal as 
in iPieri et al. (2008)| 



^""^"(i/r, AQ) 



Jm 



dM 



\ dc \ \ dOdcp \ 
Jc J J&n Ji.o.f 



dA 



[psh(M, R(Re, A, If,, 0, <p)) X P(c)x 

X KZ°(M, c, r(A, A', ^, e', 0')) X J(x, y, z\A, 0, c/>)] (9) 

where AQ is the solid angle of observation pointing in the di- 
rection of observation ij/ and defined by the angular resolution 
of the detector 0; p^h is the sub-subhaloes mass and distribution 
function inside the dwarf; J is the Jacobian determinant; R is the 
galactocentric distance, r is the radial coordinate inside the sin- 
gle sub-subhalo located at distance A from the observer along 
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Fig. 10. The astrophysical/cosmological contribution <l)cosmo to 
the y-ray flux as a function of the angular distance from the 
galaxy centre, derived for Draco using the DM mass density 
radial profiles by [Lokas et al. (2005)| (solid line). Dashed, dot- 
dashed and long-dashed lines show the contribution of all sub- 
subhaloes. The dotted line shows the contribution of the DM 
halo of Draco when only 60% of its mass is smoothly distributed 
in the halo. 



the line of sight defined by iff and contributing to the diffiise 
emission; O^";™" describes the emission from each sub-subhalo. 
As in the case of the MW, psh is normalized such that 10% of 
the Draco mass is distributed in substructures with masses in 
the mass range between 10"^ and 10"^ Moraco- When integrat- 
ing over all sub-subhaloes, we end up with 40 % of the Draco 
mass distributed in lO'^ sub-subhaloes. The results of the com- 
putation of (]>™™° using Eq|9]are shown in Fig[TO]using differ- 
ent models for the concentrations parameters, namely the BzO.ref, 
5zf,ief and Bzf,5o- described in |Pieri et al. (2008)| Superimposed 
are the values of 0™'™° for the profile by [Lokas et al. (2005)| 
with 100% and 60% of the mass of Draco smoothly distributed 
in the halo. The sum of the latter contribution and the sub- 
subhalo ones should be compared with the 100% smooth halo 
(solid line). We note that the eff'ect of adding sub-substructures 
can give an enhancement of several orders of magnitude at 
large angles, where the overall flux is though too low to be de- 
tected, even in the presence of sub-subhaloes. Yet, the effect 
in the very inner parts of the halo wiU be that of decreasing 
the expected signal, and corresponding boost factor defined as 
( ^I^^ -i- ^I^^ )(E^^ is less than 1 where the 

^ 60%smooth sub-subhaloes'^ 100%smooth 

larger detectable flux is expected, that is toward the galaxy cen- 
ter. 



3.2. Boost factor due to the presence of a black hole 

So far there are only two examples of dw arf galaxies suggested 
to host a SBH. [Maccarone et al. (2005)| discussed the possibil- 
ity that the radio source near the core of the Ursa Minor dwarf 
galaxy is a SBH. They give a mass upper limit of ~ 10"* Mq. 
[Debattista et al. (2006)| assumed that the double nucleus of the 
dE VCC 128 is a disk orbiting a SBH. They derived a SBH 
mass of ~ 10^ M©. The lack of SBHs in dwarf galaxies was 
explained by |Ferrarese et al. (2006)| They found that for galax- 
ies less massive than few 10^ M© the formation of a compact 
stellar nucleus is more likely than that of a SBH. Both stellar 



nuclei and SBHs contain a mean fraction of about 0.2% of the 
total mass of the galaxy. The same conclusion was reached by 
iWehner & Harris (2006) and Cote et al. (2006)| 

Nevertheless, a value for the SBH mass of the dSph galaxies 
studied in this paper can be inferred using the Msbh - cr rela- 
tion (see [Ferrarese & Ford (2005)| for a review). Extrapolating 
the scaling law defined by SBHs in massive galaxies to the con- 
stant cr Si 10 km s"' measured in the sample galaxies, the de- 
riv ed SBH mass is Msbh ~ 1 0^ Mq. 

[Gondolo & Silk (1999)| and |Merritt (20"04)| studied the effect 
of an adiabatically accreted SBH on a cuspy DM mass density 
profile [p(r) oc with < y < 2]. They found that the SBH 
induces a central density spike described by a power-law radial 
profile with index 



Ts 



9-2r 



4 - y 

over a region of radius 
- 0.2 rsBH = 0.2 — 



(10) 



(11) 



where tsbh is the radius of gravitational influence of a SBH with 
a mass Msbh and cr is the DM velocity dispersion at tsbh- Eq.fTOl 
reflects the condition of adiabaticity, requiring the SBH forma- 
tion time to be much larger than the dynamical timescale at a 
distance tsbh- 

Assuming the cuspy DM mass density profile [y = I) given 
by |Lokas et al. (2005)| for Draco, it results 



Psbh+dm('') 



' p(rc) r <r^ 

P('-s)(7^) r^<r<r^ 
>:exp(-^) r>r, 



(12) 



Ann 



where Mdm and r\, are the total DM mass and break radius as 
in |Lokas et al. (2005 )| respectively. A core radius - 10 kpc 
was imposed to prevent the high annihilation rate which would 
destroy the spike. On the other hand, a mass density profile with 
constant-density core (y - 0) would not allow the growth of the 
spike. Therefore, it was not considered. 

The effect of a SBH on the DM mass density profile by 
|Lokas et al. (2005)] is shown in Fig. [TT] for two extreme values 
of the Ts. For cr = 10 km s"', - 10 pc corresponds to 
Msbh - 10^ Mq and rj - 10 pc corresponds to Msbh - 10^ 
Mq. The adopted value of cr is consistent with the stellar ve- 
locity dispersion measured in the dSph galaxies we studied 
( jGiknore et al. (2007)| l. It represents an upper limit for the DM 
velocity dispersion, since the stars are a tracer population of both 
the luminous and DM components. 

The boost factor due to the presence of a SBH is 



BFsBH - 



SBH+DM 



Jgal gal,sm 



(13) 



Of course, a different parametrization of the SBH properties 
would result in different values of BFbh- In particular, if we as- 
sume = lO""* - 10"^ pc we end up with a boost factor in the 
range [1,10], while if we consider = 1 - 10 pc the boost factor 
becomes of the order of [5 x 10^, 10^]. 

A SBH mass Msbh ~ 10^ Mq results in rj = 10'^ pc, and in a 
boost factor equal to 1. On the other hand,I Gilmore et al. (2007)| 
the existence of observed small stellar clumps inside Ursa Minor 
is not compatible with the presence of a spike by which they 
would have been tidally destroyed. 



10 



L. Fieri et al.: Could the Fermi-LAT detect y-rays from dark matter annihilation in the dwarf galaxies of the Local Group? 



maximal spike 
minimal spike 

Lokas et aL 




Fig. 11. The density spike induce by a SBH on the DM mass 
density profile of Draco. The initial profile is taken from 
[Lokas et al. (2005)| (solid line). The dotted and dashed lines cor- 
respond to a final density profile with a spike radius = 10"^' 
and 10 pc, respectively. 



The MAGIC collaboration has recently published 
( [Albert et al. (2008)| an upper limit on the y-ray flux 
above 140 GeV for Draco at 10 " ph cm"^ s If we 
multiply the value of Ocosmo in the direction of the cen- 
tre of Draco by a boost factor BFbh = 5 x 10^ we end 
up with (Dcosmo Ibfbh ~ 4 X 10"* GeV^ cm ^ kpc sr. 
This means that, if such a boost factor exists, either the 
DM particle mass is below 140 GeV and its annihila- 
tion products could not be observed with MAGIC, or that 
Opp < 2.5 X 10"'^ cm'* kpc"' GeV"^ s"' sr"' which is a very 
low (though allowed) value. 



4. Conclusions 

The Fermi-LAT telescope was launched in June 2008 and is tak- 
ing data on y-rays in the energy range between about 20 MeV 
and 300 GeV. Its all-sky survey operation mode will allow an 
unprecedented precise study of the y-ray sky, so that many DM 
models will be tested. The dSph galaxies of the Local Group will 
be primary targets for DM analysis with the Fermi-LAT, because 
of the low astrophysical background expected in their direction. 
Therefore, we studied the detection limit of the y-ray flux from 
DM annihilation in four of the nearest dSph galaxies, namely 
Draco, Ursa Minor, Carina, and Sextans. 

State-of-art DM density profiles were available for these 
galaxies and we computed the expected y-ray flux from DM an- 
nihilation for different particle physics parameters. We varied the 
DM particle mass as well as the annihilation cross-section and 
branching ratios. We found that the presence of NFW-like cusp 
or constant density core in the DM mass density profile does not 
produce any relevant effect in the y-ray flux due to a combina- 
tion of the geometrical acceptance of the Fermi-LAT detector, 
which is not able to resolve the very inner shells of the studied 
galaxies, and the distance of the sample galaxies. 

In the case of Draco and Ursa Minor, we found that they 
would shine above the Galactic smooth halo for all but the small- 
est angles (~ 2 degrees) above the Galactic Center Yet, the upper 
values of the predicted flux were found to be about two orders of 
magnitude below the Fermi-LAT detection threshold as derived 



in [Baltz et al. (2008)] Such values were computed for the most 
optimistic particle physics scenario of a 40 GeV particle with 
Canni' = 3 X 10"^^ cm"' s"' annihilating into bb. We have shown 
how the effect of the boost factor due the presence of a popu- 
lation of DM clumps inside the dSph galaxies, though possibly 
dramatic (of the order of 10^ when integrated over the whole 
galaxy), had to be rescaled for the limit on the EGB measured by 
EGRET. The overall maximal effect was reduced to a factor 70. 
The reader should keep in mind though that the calculation was 
made for a toy-model where the subhalo population of the dwarf 
galaxies is described by the same model used for the MW. In any 
case, since the closest dwarfs are not pointlike for the Fermi-LAT 
angular acceptance Af2, the factor to be taken into consideration 
is the effect of the sub-subhaloes inside Af2, which resulted in a 
decrease of the expected flux due to a redistribution of the DM 
inside the halo. The presence of a central SBH in agreement with 
the MsBH - cr relation extrapolated to the observed low cr val- 
ues resulted in a negligible boost factor More extreme models 
would result in a much higher boost factor, though they are not 
theoretically supported. 

Contrarily to previous papers which addressed the presence 
of subhaloes or of SMBHs to boost the signal, we have demon- 
strated that the boost factor must be searched for in some exotic 
extension or modification of the particle physics sector Unless 
such a scenario happens, we conclude that the annihilation of 
DM inside the local dwarfs is unlikely to be detected with the 
Fermi-LAT. 

As a further development, it will be interesting to repeat the 
present analysis of boost factors for the recently catalogued ul- 



tra faint dwarfs of the Local Group. Strigari et al. (2008) have 



computed the expected y-ray flux from those sources, deriving 
the halo parameters from kinematical data. The inclusion of such 
galaxies in our study will improve the sensitivity of a joint multi- 
centred likelihood analysis. 

It is worth noticing that,since the DM spectra would be the 
same for all the galaxies, such an analysis could be performed in 
order to maximise the detection efficiency and to allow portions 
of the particle physics phase-space to be explored. 
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